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Abstract 

In the Reflection Positivity theory and its application to statistical mechanical systems, certain 
matrix inequalities play a central role. The Dyson-Lieb-Simon [1] and Kennedy-Lieb-Shastry [2] in¬ 
equalities constitute prominent examples. In this paper we extend the KLS inequality to the case 
where matrices are replaced by certain operators. As an application, we prove the occurrence of the 
long range order in the ground state of two-dimensional quantum rotors. 
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1 Introduction 

The Reflection Positivity notion has appeared in Quantum Field Theory in seventies of the last century 
[3]. Few years later, it has been applied to investigation of phase transitions in both classical |1] and 
quantum [T] lattice spin systems. The Reflection Positivity turned out to be a very useful tool, giving 
the first rigorous proofs of existence of phase transitions in systems with continuous symmetry group. 

The cornerstone of Reflection Positivity for quantum spin systems is the matrix inequality due 
to Dyson, Lieb and Simon (Lemma 4.1 in [T]). Using this Lemma, authors proved the existence of 
orderings in the XY as well as Heisenberg models in c? > 3 and for sufflciently small temperature. 
Later on, this method has been extended to certain class of infinite dimensional operators. This way, 
the existence of Long-Range Order has been proved for d > 3 in the system of quantum interacting 
rotors |S]. 

Another direction of development of Reflection Positivity techniques was an examination of ground 
states of quantum spin systems and orderings therein. It turned out that one can take certain zero- 
temperature limit in the framework of the DLS method. This way, the appearance of Long Range 
Order has been proved in XY and Heisenberg models in d = 2 0 , 17 ], 0 . Later on, it turned out that 
such a proof can be done directly in the ground state, with the use of another matrix inequality, due 
to Kennedy, Lieb and Shastry (KLS) 0. This inequality was further generalized by Schupp 0. 

It would be tempting to extend this inequality to infinite-dimensional version, i.e. for certain class 
of operators. However, to our best knowledge, the operator version of the KLS and Schupp(KLSS) 
inequalities, suitable for applications to ground states of quantum interacting rotors has not been 
developed. 

This opportunity inspired us to attempts to prove an operator analog of the KLSS inequalities. It 
turned out to be possible, and this is one of two main results of our paper: extension of the KLSS 
matrix inequalities to certain class of infinite-dimensional operators. The second group of results which 
seems to be new are some applications. 
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The outline of the paper is as follows. In the Sec. [^we formulate the operator version of the KLSS 
inequalities. The application of this inequality is described in the Sec. it is the proof of the ordering 
in ground state of d > 2 rotors (alternative proof to that presented in [10] )■ The Sec. [^contains 
summary, conclusions and description of some open problems. 


2 KLS inequality and its extension for operators 


2.1 Kennedy, Lieb, Shastry and Schupp matrix inequalities. 

For convenience of the reader, and to show the idea of a proof without operator-theoretic details, we 
present firstly the matrix version of KLSS inequality. 

Theorem 2.1 |5| Let c,A^B be n x n complex matrices, |c| := Vc*c and |c*| := '/cc* the moduli of 
c and c* respectively. Then 

\Tyc* BcA*\ < i[Tr (|c| A |c| A*)Tr {\c*\B\c*\B*)]. (1) 

Sketch of the proof: At first let us note that by the polar decomposition theorem c is of the form 
c = u |c|, where u is a partial isometry. Since u* u |c| = |c| and u |c| u* is a positive matrix, the polar 
decomposition of c* is of the form: 


c = \c\u =u u\c\u =u c 


Taking adjoint we get c = u |c| = |c*| u. Therefore 

w\/R = \/\^\u 


( 2 ) 

(3) 


according to functional calculus of positive hermitian matrices and 

c=\/|c1w\/M- (4) 


Now, let P and Q be matrices introduced by formulae 

P = u* ^\B* ^\u, Q=^\A*yM- (5) 

Let us remind that the trace functional defines the scalar product on the space of square matrices: 
{A\B) := Tr A*B. Therefore 

|TrP*Q| < y(Tr P*P) (Tr Q*Q) < ^[Tr + Tr Q*Q] (6) 

due to the Schwarz inequality followed by inequality between geometric mean and arithmetic one. 
Now using formula @ and (@ one can easily verify that 

P*Q = u* y/\(^\ B y/\^\ u\/\c\A* ^/\c\ = '/\c\u* Be A* y/\c\ 

and due to 0 

Tr P* Q = Tr = Tr {\c\u* Be A*) = Tr {c* BcA*). 

The module of it coincides with the left hand side of the inequality 0. 

In the similar manner we compute the right hand side of 

P*P = u* v1^ B uu* B* v1^ u = u* B\c*\B* ^/\P\ u. 


Therefore 

TrP*P = Tr (^u* ^/\P\B\c*\B* ^/\P\u^ =Tv {\c*\B\c*\B*). (7) 

By the similar reasoning we get 

Tr Q* Q = Tr ( A\c\ A* = Tr (|c| A |c| A*). (8) 

Combining ^ and ^ we obtain the right hand side of inequality Q. 
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Theorem 2.2 m The KLS inequality 0 holds also for rectangular matriees, i.e. c is n x m matrix 
and matrieies A and B are m x m and n x n respeetively. 

Proof: It is almost a repetition of the proof for KLS inequality and refers to the modified polar 
decomposition of the operator c acting between m-dimensional and n-dimensional spaces. In that 
situation |c|, |c*| are positive matrices of dimensions mx m and n x n, respectively. ■ 


2.2 Operator version of the KLSS inequality 

The main goal of this subsection is to prove a generalized version of the KLSS inequality: 

Theorem 2.3 Let C and TZ be separable Hilbert spaces, A € B{C), B € B(7Z) bounded operators 
acting on C and TZ respectively. Let c : C TZ be a Hilbert-Schmidt operator and |c| := V c*c, and 
|c*| := y/cc* be the corresponding moduli. Then 

1. |c| and |c*| are hermitean Hilbert-Schmidt operators acting on C and TZ respectively; 

2. c*BcA*, |c|A|c|^* are trace-class operators on C and |c*|i3|c*|i?* is a trace-class operator on 

TZ; 

3. the following inequality holds 

|Tr c* BcA*\ < ^ [Tr (|c| A |c| A*) + Tr {\c^\B\c*\B*)] (9) 

Remark. In a finite-dimensional case, i.e. C = — TZ the above inequality reduces to the matricial 

KLS inequality [2] and in more general finite dimensional situation dim C ^ dim TZ we obtain the result 
of Schupp 12. 

Proof: To prove our result we shall use some properties of Schatten ideals m, CH], m and now we 
shall recall necessary results of the theory. 

Let TLhe a separable Hilbert space, CB{TL) - the set of compact operators on TL. For a real number 
p>l the p-Schatten ideal is the set 

LP{H) := {a G CB{TL) : Tr (|a|P) < oo}. 

Remark. Sometimes LP{'H) is denoted as ffp jTJj, however, the actual notation corresponding to 
noncommutative L^-spaces seems more natural to us. 


For a G LP{fH) let us define: 

||a||p := (Tr|ar)-. 

Then it is known that 

i) {LP{TL), II • lip) is a Banach space, LP(H) is a two-sided ideal in BifH), i.e. for any a G LP{H) and 
A, B G B{H) the operator AaB G LP{Ti) and moreover 


||Hai?||p< ||H||||R||||a 


a) A p,q,r > 1 are such numbers that: —I— = - and a G LP{TL), b G L’^LH), then ab G L'^iTL). 

p q r 

Hi) in particular, ifp and q satisfy - -f- - = 1 then for a G LP{fH), b G L‘>(H) products ab, ba G L^ifH), 

P q 

Tr ab = Tr ba and 

|Tr a6| < ||a||p ||6||,. 

Remark. The space L^{'H) is the space of trace-class operators on TL and Lf{TL) is the space of 
Hilbert-Schmidt class. Clearly Lf{TL) equipped with the sesquilinear form 

L’^i'H) X L^{Ti) 9 {a,b) i— ^ (a|6) := Tr a*b G C 


is a Hilbert space. 
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In what follows we shall also need Hilbert-Schmidt operators in more general settings, namely the 
operators from one Hilbert space to another. 

Definition. An operator c : £ —> 7?. is a Hilbert-Schmidt one, if for some orthonormal basis {a} 
in C, the sum 

I|c|l 2 := '^{ca\ca) = ^ \\ca\\^ 

OL Ol 

is finite. 

The set of such operators will be denoted by L^{C, TZ). Clearly for c S TV) and any orthonormal 

basis {/?} in TZ, we have 

J2{ca\ca) = ^(a|c*ca) = ^ |(/3|ca)P 

CK a q ;,/3 

= ^ |(c*/3|a)p = ^ ||c*/?|p = Y.mcc*P). (10) 

a,/3 0 0 

Therefore c* S L'^{TZ,C) and the finiteness condition does not depend on the particular choice of an 
orthonormal basis {a}. In particular c*c G L^(£), cc* € L^{TZ) i.e. they are trace-class operators 
acting on C and TZ respectively. Moreover 

||c||2 = ||c*c||i = Tr£(c*c) = Tr,^(cc*) = ||cc*||i = ||c*||2. (11) 

Let us note that for a, b G L‘^{C,TZ) we have a* b G L^{C), ba* G L^{TZ). Moreover one can easily 
check that for any Hilbert spaces C', C,TZ,TZ' a modified ideal property holds: 

/ cGL^{C,TZ) \ 

and ^ {BcAgL^{C',TZ')) . (12) 

V AgB{C',L), BgB{TZ,TZ') j 

As before the space L‘^{C,TZ) forms a Hilbert space equipped with the scalar product 

(a|6) := Trcia* b) = Trn{ba*) = {b*\a*). (13) 

The last equality can be verified by the similar calculation as above. In what follows to simplify 
notation the corresponding indices C or TZ will be omitted. As the result, by the Schwarz inequality 
followed by mean arithmetic-geometric inequality, we obtain 

Corollary 2.4 For arbitrary a,b G Lf{C,TZ) we have 

|Tr a*6| < i [Tr a*a -k Tr 6*5] (14) 

Now we are ready to prove our result. Assume that c G Lf{C,TZ). Therefore c* c and cc* are 
trace-class i.e. jcj G and |c*| G li^iTZ) and this proves the first part of the theorem. The second 

part easily follows from ( |I^ . 

To prove the inequality let us note that by the polar decomposition theorem c is of the form 
c = u jcj for the unique partial isometry u G B{C,TZ), u : |c|(£) -G c{C) such that u* u and uu* are 
projections on the initial and final domain respectively. Now, by uniqueness of the polar decompostion 
and functional calculus of bounded, self-adjoint operators we obtain (in the same way as for matricies): 


U'Jic] = ^/\^\u. (15) 

Therefore _ 

C= ^/\c^\uy/\c\. (16) 

Clearly \/[c[ G £^(£) and -^/jc* | G C^{TZ) and this observation enables us to follow the proof given 
for matrices in 12 ], i- As in ([^ we define operators P,Q G B{C): 

P = u* vl^w, Q = v^A*v1^. (17) 

Remembering that L^-spaces are ideals and using pr ope rty ii) of Schatten ideals (f or p = g = 4) we 
see that Q G T^(£) and y/\c* \ B* yj\c* \ G L'^{TZ); by (12) we have P G ■ Now (14) reads 

\Tr P*Q\<^[Tt P*P + Tr Q*Q] (18) 
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To compute the left hand side of the above expression let us notice that u* by 

formula (15l. Using this fact and (16 1 we get 

P*Q = u* y/\c^\B ^/\c^\u ^\A* ^\ = ^\u* B cA* ^\. 


The operator u* Be A* belongs to L^{C) due to the modified ideal property (12 Since \/\c\ G L'^{C) 
the operator ^/\c\ u* B c A* G (£) by property ii) of Schatten ideals. Now using property in) in the 


case P = I and g = 4 we have 


Tr = Tr (^^\u* BcA* = Tr {\c\u* BcA*) = Tr {c* BcA*) 

due to ([^ . The module of it coincides with the left hand side of the inequality (|^ . 

In the similar manner we compute the right hand side of (18 1 . 


P*P = u* B \/i^ u u* a/H B* a/H u = u* B\c*\B* a/Ic^ u. 


Therefore 


Tr P* P = Tr 


(u* ^\B\c*\B* = 11 ^^ 

^ = Tt (y^\B* ^\uu* 


= Tr [^\B*\c*\B^\) 


by (11). Now a/|c*| B* G L^iJZ) and |c*|P € La{C) and using property in) again we obtain 

Tt P* P = Tr {\c*\B\c*\B*). (19) 

By the similar reasoning we get 

Tr Q* Q = Tr A\c\ A* ^\'^ = Tr (|c| A |c| A*). (20) 

Combining (19) and (20) we have the right hand side of inequality ([^. The proof is done. 


2.3 Main inequality and expectation values 

In this subsection the main inequality will be expressed in terms of expectation values of operators 
acting on £ 0 7Z, where £ and 7Z are separable, infinite dimensional Hilbert spaces. In that form it 
will be used in following sections. 

Let r := {V' 7 } denote a fixed orthonormal basis in £. It defines a linear map: 


f : L^(£, P) 9 c !->■ f (c) = ^ V't C) c-ipj G £®TZ 


Basic properties of this map are described by 

Lemma 2.5 The map T is unitary, moreover for B G B(TZ) and a,cG Lf(£,TZ) we have: 

f(Pc) = {I®B)f{c) 


(f(a) I (Jo P)f(c)) = (f(a) I f(Pc)) = Try(a*Pc) = Tr 7 ^(Pca*) = Tr 7 ?,(ca*P) (21) 

Proof: It is straightforward to verify that T is an isometry i.e. for Ci,C 2 G Lf{£,7V) we have 
(ci)|f( 02 )^ = Tr£(c)(c 2 ). The formula for f(Pc) is clear. To show unitarity, notice that any 
/i € £ O P is of the form h = ® for the unique family of vectors (r^) in P. Let us define a 

linear map ri(/i) : £ —)■ P by Ti(h)4’--i '■= r^. Then ri(J) G £^(£,P) and simple calculation shows 
that rri(/i) = h and rir(c) = c, so r is unitary and Ti = T*. The formula (21) is also clear; the last 
two equalities follow from properties of trace: Tr£(a*c) = Tr 7 ^(ca*). ■ 
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The basis T defines also an antiunitary involution J-p '■ C ^ C: 


Jr := kjtljy , fc-y G C ; Jp = Jr ; Jy = I. 

Note that for A G L^{C) we have: 

Tr£(JrAJr) = ^ ('0-y | JyAJyljjj) = ^ {AJytpy \ JpV't) = ^ I V't) = 

= Trc(^*) 

In the same way for a basis fl = {^(^} in 7?. we have the mapping f2 


fl : (7?., C) 9 J i-G fi(J) = ^ d(j)aj ® 4>ui G C®TZ] 

with the corresponding antiunitary involution Jq : Jq :=^ L4>uj, and 

Lemma 2.6 T/ie map 17 is unitary, moreover for A G B{C), B G L^{TZ) and b,d,G L‘^{TZ,C): 

n{Ad) = (A (g) i)n{d) 

in{d) I {A (g) I)n{b)) = (n{d) \ n{Ab)^ = Tin{d*Ab) = TrdAbd*) = Trcibd*A) 

TTniJnBJn)=TTn{B*) 


(22) 


(23) 


(24) 


The choice of bases in C and TZ gives us all of these objects and the following straightforward 
lemma describes relations between both structures: 

Lemma 2.7 Let c G L‘^{L,TVj and d G L^{TZ,C). Then the following equalities hold: 

f(c) = i7(JrcVo), f(JaJVr) =i7(J) (25) 

(f2(J) I f(c)) = (n{d) I f7(Jrc* Jn) 

= (f(JaJVr)lf(c) 

■ 

Clearly, the choice of bases T, LI is equivalent to the choice of a basis T and a unitary operator 
U : C ^ TZ (the equality U{ij}.y) = (l>j is a definition of operator C7 or a basis LI = {(fy})- It is sometimes 
more convenient to use pair (T, U) instead of (T,!!). In the lemma below we collect formulae we will 
use: 

Lemma 2.8 For c G L‘^{C,TZ), A G B{C), B G B{TZ) the following equalities hold: 


Jn 

= UJrU* 

(27) 

(f(c)|(J®i?)f(c)) 

= (f(|c*|t7)|(J®i?)f(|c*|C)) 

(28) 

(f(c)|(^®/)f(c)) 

= (f(C/|c|)|(Jl®/)f(t/|c|)) 

(29) 

(f(c)|(^®/)f(c)) 

= (f(C/|c|) 1 (J® 7/Jr^*JrC/*)f(C/|c|)) 

(30) 

(f(c)|(J®i?)f(c)) 

= (f(|c*|t7) 1 {JrU*B*UJr 0 /)f(|c*|C/)) 

(31) 

(f(c)|(Jl®i?)f(c)) 

= Treic* BcJyA* Jr) 

(32) 


= Tm{d*Jrc*Jn) = 
= Tr£(JrJJnc) 


Proof: For c G L'^{C,TZ) recall that |c| := y/c*c G Lp‘{C) and |c*| = \/cc* G L‘^{TZ); notice also that 
U\c\ G L^{C,TZ) and |c*|I7 G L^{C,TZ). 

The formula (27) is straightforward, let us prove (28). For B G B{TZ), compute: 

(f(c) I (7® B)f(c)) = Tide*Be) = TTnicc*B) = Trni\c*\^B) = Tr 7 ^(|c*|i?|c*|) = 

= TrciU*\c*\B\c*\U) = Ti c{{\c*\U)* B\c*\U) = 

= (f(|c*|C7)|(/®i3)f(|c*|C/)) 
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For the next formula, let A € B(C) and using lemma 2.7 we compute: 


(f (c) I (A 0 /)f (c)) = (fi(JrcVa) | (A ® /)fi(Jrc* Jn)) = TrK((JrcVn)*AJrcVa) = 
= Ti'7?,('^acJrA( Jpc* Jo)) = Tr£( Jpc* Jo JocJrA) = 

= Tr£(Jr|c|VrA) 


Writing the formula above for U\c\ instead of c and noting that |C/|c|p = (C/|c|)*C/|c| = |cp we obtain 
the equality (29l. 


We prove the equality (30): 


(f(c) I (A(g) /)f(c)^ = Tr£(Jr|cpJrA) = Tr£((Jr|c| Jr)Jr|c| JrA) = Tr£(Jr|c| JrAJr|c| Jr) = 
= Tr£(|c| JrAVrIcI) = Tr£((C/|c|)*C/Jr AVrC/* (C/|c|)) = 

= (f (C/|c|) I (7 (g) C/JrAVrC/*)f (7/|c|)) 


and the formula (31l: 


(f(c) I (7® S)f(c)) = Tr 7 ^(|c*|S|c*|) = Tr7^(Jo|c*|B*|c*| Jo) = 

= Tr,^(Jo|c*|7/Jr JrC/*B*C/Jr JrC/*|c*| Jo) = 

= TrK((JrC/*|c*|Jo)*(JrC/*B*7/Jr)(JrC/1c*|Jo)) = 

= (^h{JTU*\c*\Jn)\{JTU* B*U Jt ® m{JTU*\c*\Ja)) = 
= (f (|c*|C/) I {JrU*B*UJr 0 7)f (|c*|C/)) 


Finally, we prove (32l: 


(f(c)|(A®i?)f(c)) 


((A* 0 7)f(c) I (7 0 B)f(c)) = ((A* 0 7)f)(Jrc*Jo) | f(Sc)) = 
(fi(A* Jrc* Jo) I f (Sc)) = Tr7j( JocJrAJrc*S* Jo) = 

TrK(BcJrA* Jrc*) = Tr^) Jr A* JrcVo Jo^c) = 
Tr£(JrA*Jrc*Sc) = Tr£(c*ScJr A* Jr), 


where we have used (261. 


Now we want to express the inequality @i.e. 

2|Tr£(c*BcA*)| < Tr£(|c|A|c|A*) + Tr 7 ^(|c*|B|c*|S*) 
in terms of F, Jr, U and the scalar product in £ ® 7?.. 


(33) 


In the formula (32) we put B := UJrAJrU* and U\c\ instead of c and get: 


(f(C/|c|) I (A® (7/JrAJr[/*))f(C/|c|)) = Tr£(|c|C/*(t/JrAJr7/*)(C/|c|)JrA*Jr) 

= Tr£(|c| Jr AJr |c| JrA* Jr) 


Now put into (32) A 


U* JnBJall and |c*|77 instead of c and obtain: 


(r{\c*\U)\iU*JnBJnU®B)ri\c*\U)^ = Ttc{U*\c*\B\c*\UMU* JnBJnU)* Jr) = 

= TYc{U*\c*\B\c*\{UJrU*)JnB*JnUJr) = Trc{U*\c*\B\c*\{JnJa)B* Jn{JnU)) = 

= Ttc{U*\c*\B\c*\B*U) = 

= TTn{\c*\B\c*\B*), 


where (27) i.e. Jo = UJrU* was used. 

Finally, writing (33) with J^AJr instead of A, using (32) and two equalities above we obtain: 
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Proposition 2.9 Let A G B{C), B G B{TZ) and c G L‘^{C,TZ). Then 


+ 


(f(c)|(^®B)f(c))| < (r{U\c\)\{A<^UJrAJrU*)r{U\c\)^^ 

+ (f(|c*|C/) I {U*JnBJnU<g)B)r{\c*\U)) 


(34) 


We will also need formulae similar to the ones in Lemma [2l8| for some unbounded operators, they are 
proven in the following proposition. 

Proposition 2.10 Let T and S be self-adjoint operators with purely point spectrum acting on C and 
TZ respectively. Assume that bases {if'y} o-n-d {0^} consist of eigenvectors of T and S: Tijjj = tryipry 
and S(j).y = Sj(j)j; assume moreover that r(c) G D({T ® I)^) and n{d) G D({I ® 5')^). Then: 

(f(c) I (T 0 /)f(c)) = (f(C/|c|) I (T 0 7)f(C/|c|)) = (f(C/|c|) | (J ® UTU*){f{U\c\)) (35) 

(fi(d) I (7 (g) s)h{d)) = (^(c/* Ml) I (7 (g) s)h{u 


(n{U*\d\)\{U*SU:^T)n{U*\d\)'^ (36) 


Proof: We will prove (35); equalities in (36) can be proven in a similar manner. 

Since r(c) is in the domain of (T (g) 7)^ we have: 

(f(c) I (T 0 7)2f(c)) = ® I (T ® Ifr{c)) = ® c ^7 I f (c)) = 

7 7 

= i7nicV'7ir 


(37) 


This equality means that the series is convergent. Because r(g7 is closed and G 

D{T (g) 7) it implies that: 


(T (g) 7) f (c) = ^ tjlf-y (g) Clf-J 


(38) 


Since ||C/|c|'(/' 7 lP = ||c^.y|p the formula (37) implies also convergence of the series ® Ulclip. 

and the equality 


(T®7)f(C/|c|)=^ ^'07 (g) U\c\lp-y 


(39) 


By the lemma 


2.7 


r(?7|c|) = n{ Jr\c\U* Jq). Since Jr|c|t /*(/>7 ® = Jr|c |'!/'7 ® 4>^ G D{I ® UTU*) 

the formula (37) means convergence of the series ’^.ytjJrlclfjj ® (j).y and the equality 


(7(g) C/TC/*)f(C/|c|) = ^^(Jr|cM7) (g ^7 


(40) 


follows. Now combining (38), (39), (40) and r(t/|c|) = LI{Jy\c\U*Jq) we obtain (35). 


3 Ground state ordering in the system of 2-dimensional rotors 

3.1 Description of the system 

Denote by A the finite subset of the simple cubic lattice in d dimensions: A C We assume that A 
is a (discrete) hypercube and that the number of sites along every edge is even] let us fix 2N to be the 
length of the hypercube edge: 

A := {x € : —N 1 < Xi < N, i = 1,... ,d} (41) 

With every site x G A we associate a real variable (px G [0, 27r[. In physical terms, it describes the 
(angular) position of the rotor at the site x. Equivalently, the position of the rotor at the site x can 
be described as a unit vector Sx G S^, i.e. one dimensional torus T: 

Sx = [Sx:4] = [cos(px,sin(px]- 
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The total spin is S = Sx- The Hilbert space "Hx of states on a given site x is the space of square 

xGA 

integrable periodic functions, i.e. "Hx = T^(T). The Hilbert space associated to the whole system is 
the space of square integrable functions on |A| - dimensional torus: 

•Ha = = 0 Hx 

xGA 


The operator T of total kinetic energy of the system of rotors is proportional to the laplacian A: 


T 


1 

AT 


E 

xGA 




(42) 


where / > 0 is the moment of inertia of rotor (we assume that all rotors have equal moments of inertia). 
The system of interacting rotors is defined by the Hamiltonian H = T + y, where V is an interaction 
energy between rotors; it is an operator of multiplication by a smooth function V. We shall consider 
the Hamiltonian: 

H = T + yo, Vo:=-jJ2cosiip^-ipy) (43) 

(xy) 


In this formula (xy) means that x and y are the nearest neighbours. By this we mean that all but one 
coordinates ofn. and y are equal and the ones, say Xi and yi, that differ satisfy \xi—yi\ = 1 mod 2(7V—1). 
Thus a site x laying on the hyperplane defined by xi = N has some of its nearest neighbours on a 
hyperplanes defined by Xi = —N + 1. 

J is the coupling constant: J > 0 corresponds to ferromagnetic coupling between rotors and J < 0 
to the antiferromagnetic one. In what follows we shall restrict ourselves to the ferromagnetic case, as 
only in this situation the Reflection Positivity arguments may be applied. 

Our Hamiltonian (431 is an elliptic second order differential operator on a |A|-dimensional torus 
(compact manifold). It is a special case of the more general situation: 


Theorem 3.1 Let (M,g) be a compact, oriented, riemannian smooth manifold without boundary; 
L : C°^{M) —)■ C°°{M) formally selfadjoint, linear, elliptic, PDO of order k > 0. Then: 

1. L extends uniquely to L : H^{M) —> Lf{M); is k-th Sobolev space) 

2. L as an operator on Lp'{M) (with the domain H^{M)) is selfadjoint; 

3. The spectrum of L consists of isolated eigenvalues of finite multiplicity; 

4 . Eigenvectors of L are smooth functions. 


Let us observe that since T is positive and Vq is a continuous function, the hamiltonian H is bounded 
from below. 

In the following we need in an essential way the uniqueness of the ground state of the Hamiltonian. 
This is a consequence of positivity improving property of the semigroup exp(—tA). For convenience of 
the reader, we recall briefly main definitions and results (we refer to Chapt. XHI of [16] for a detailed 
presentation). 

A non zero function 4* G L‘^{M) is positive iff 4'(a;) > 0; it is strictly positive if 4'(a;) > 0 (both 
inequalities should be understood in almost everywhere sense). 

A bounded operator A is: 

- positivity preserving if A'^ is positive for positive 'k; 

- positivity improving if A4' is strictly positive for positive 4'; equivalent condition is that (4* | A4>) > 0 
for positive 4* and 4>. 


The following is, simplified for our needs, Thm XHI.44 from m- 


Proposition 3.2 Let H be a self adjoint, bounded from below operator on L‘^(M). Assume the spec¬ 
trum of H consists of isolated eigenvalues of finite multiplicity. If for every t > 0 the operator exp{—tH) 
is positivity improving then the ground state of H is unique (and strictly positive). 


In our situation M = and it is known that for kinetic energy operator T given by the formula 


(421 operators exp(—tT), t > 0 are integral operators; for $, 4^ G L^(TI^I) : 


($ I exp(-tT)4') 


J dxdy^{x)K{t,x,y)'I>{y). 
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The function K{t, x,y), t > 0,x,y G - the heat kernel for |A| - torus, is explicitely known; it is a 
strictly positive function, therefore exp(—tT) is positivity improving. 

The same property for H = T + V holds due to the Trotter product formula (see e.g Chapt. X of 

my- 

'' r 1 ^ 

exp(—t(T + y))^'= lim exp(—tT/n) exp(—ty/n) dr 

n—>-oo L 

In fact, for any $, dr S Ha : 

(d) I exp(—tiJ)dr) = lim (d)| exp(—fT/n) exp(—ty/n) 

n—^oo \ L 


d- 


Let C := supy. Since exp(—tT) is positivity improving, for strictly positive numbers: a,(3,"f,d > 0 
and positive functions $ and dr we have: 


d) I exp(—aT) exp(—/3y) exp(— 7 T) exp(—(5y)dr^ > e ^exp(—Q;T)d’| exp(— 7 T) exp(—(5y)dr^ 

> e-y+'5)C I exp(-(a + 7 )T)d'). 

In particular ^d> | exp(—tT/n) exp(—<y/n) 


> 


drj > (d> I exp(-tT)d') > 0. 

Therefore the inequality is preserved in the limit: (d? | exp(—tiL)dr) > 0 and the uniqueness of the 
ground state of H follows. 


3.2 Criteria of ordering 

The simplest definition of the order parameter would be an average of the total spin. However, 
this definition is of little use for the zero field (i.e. as a measure of the spontaneous magnetization) 
as it is zero due to symmetry. The more physical definition is a zero-field limit of magnetization: 
Ad = limM(/i) {h denotes magnetic field) but it is difficult to deal with. More easy to handle is the 

/ i —>-0 

average of the square of spin. It follows that if the average of the square of spin is different from zero, 
then the zero-field magnetization is non-zero, too (Griffiths theorem - see [T]). So, we take the average 
(S^) as a measure of order parameter. All averages considered in this paper are taken over the ground 
state. Following this idea, we will prove that if / and J are sufficiently large the ground state of system 


of interacting rotors, described by the potential (431, exhibits Long-Range Order (LRO): 


Theorem 3.3 Assume that I and J satisfy the inequality 


'/IJ > 


dk 


(27r) 


^/W) 


= X, 


(44) 


where the function £ : 


is defined by: 


d 

£{k) = d— coski. (45) 

i=l 

Then there exists C > 0 such that, for sufficiently large |A|.' 



Remarks. 


1. The estimation of the form similar to the one in (46 1 appeared in P-based proofs of LRO in other 
classes of models including classical and quantum spin systems n, g], 0. 


2 . The integral Ii is divergent, so above theorem does not prove existence of ground-state LRO 
in d = 1 but scaling and field-theoretical arguments [ 12 ] show (non-rigorously) that there is no 
ordering in this case. 


3. For d> 2 the integral Id is finite , so there is LRO in the ground-state. Numerical values of Id 
for physical dimensions are: I 2 « 0.909173; I 3 « 0.643954. It indicates the tendency to ordering 
increases with the growth of a dimension. 
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4. For positive temperatures, there is no ordering In dimensions 1 and 2. This statement, analogous 
to the famous Mermin-Wagner theorem for spin systems [13], has been proven in |3]. 


5. It is known that in d > 3, the LRO is present in sufficiently low temperatures, and so in the 
ground state, too |S|. Therefore the Theorem 3.3 is most interesting in d = 2. Such a result, 
using another RP arguments, has been proven in cni. In different approach, by scaling and 
field-theoretical arguments, it was obtained (non-rigorously) in |12| . 


3.3 Estimations for basic functions 

To prove the theorem |3.3| it will be convenient to work with Fourier-transformed spins. Let us define: 


a zk x 
C- 


{a = x,y) 


(47) 


where k takes value in the first Brillouin zone, i.e. kj £ i ■ • ■ > > tt} for j = 1,..., d. 

Let us remark that due to the symmetry of our system we have: 

Lemma 3.4 For the system of rotors described by the hamiltonian operators sJJ (a = x,y) satisfy 
identities: 

{sZ)=0 and = (48) 

Proof: Let Rg be the rotation (in all variables) by 9: 


{Ref){<pi ,..., i^|A|) := f{<Pi -0,..., ip\A\ - 0) 

It is clear that Rg is unitary. Since the hamiltonian ( |43| ) commutes with Rg and the ground state is 
unique we have Rgtpo = Xifo for a complex number A with |A| = 1. Therefore 


(Al) = (V'o I AlV'o) = {Rei’o I ARgfjo) = {R*gARg) 


It is easy to check that for any k: Rt^s^ = —s^Rtt and i?,r/2Sk = Now equalities (48) are 

clear. ■ 

Let us denote by g\^ the two-point correlation function in the momentum representation: 


dk = (sl(s£)* 


(49) 


Clearly g\^ > 0. Since is normal g\^ = {{sff)*s^). Remember also that depends on parameters 
/ and J. The gk function and its estimation will play a crucial role in the proof of existence of 
spontaneous magnetization. 


With this notation the inequality (46) can be rewritten as 


> C > 0. 


(50) 


The strategy of the proof of Theorem 3.3 can be described as follows (general ideas are similar to ones 
in ). First (and rather easy) step is the equality 




(51) 


Indeed, for the groud state ipo of id: 

Tl^k = I s^V’o) = I cosq^ye^'^y-ifo) = I 

k k k.x,y ' ' X 

where we have used the obvious formula: = |A|(5(x — y). In the same way we obtain: 

k 

X)k(®kV'o I SkV^o) = d’xV'o I V'o)- Adding these two equalities we get: 

Y.{su§^^r) + {si{sir) = \A\ 

k 


11 









and the equality (511 follows due to the second formula of (481 

It turns out, that for k yf 0, and this is the place where the RP arguments and the main inequality 

(52) 


(34) is used, the function can be estimated from above by an integrable (for d > 1) function: 

1 


9k < 




k yf 0. 


where the function £(k) was defined in (45l. 
With (511 and (521 in hand we can write 


1 1 

W\S°-o-2^9k> 


Now the sum 


inequality 


|A| 

(27r)'^ 1 




2 VIJ 


V |A| yyy 


E 

k=jiO 




converges as A ^ to the integral f 

J f— tt-ttI' 


dk 




. Therefore the 


implies that there exists C > 0 such that for sufficiently large A we have 




> c 


i.e the estimate (46). This way to complete the proof of the theorem 3.3 it remains to prove the estimate 
(52). The existence of such an estimate seems to be a quite general phenomenon, but at present we 
can prove it only using RP techniques. 

The inequality (461 can be viewed as the appearance of the macroscopic occupation of the k = 0 
mode. In the other words, it is an indication that in the thermodynamic limit, the g\^ function possess 
non-zero 5 function contribution at k = 0. 

We will get the inequality (52l by relating g\^ to other functions, in particular susceptibility (57l, 
for which we will get an estimate by RP techniques. 


Let {'4’n)n>o be an orthonormal basis consisting of eigenvectors of the Hamiltonian (43) with Hifn = 
En'f’n and 'j/jq be the ground state. We have: 


5k = {{sl)*sl) = (s^V'O I sllfo) = y {sltfo I V'n) (V'n | = 

n>0 

= E I I '^kV’o) = E K^kV’O I ■ 


(53) 


n>0 n>0 

There is no term with n = 0 due to the first formula of ( |48[ ). Since is normal be writtes as: 

(KSkV^OlVnJI + IRSkj lyo I 

n>0 

Next, we define: 


5k = {slislY) = y |((s(J)>0 I ^n)l^ = ^ y [l(s£V'0 I Yn)? + |((s£)*V'0 I Yn)? 

n>0 


(54) 


'Dk :=^([[5£,i^],(OT 


(55) 


Since and {sff)*4’o are in the domain of H this definition is correct Short calculation shows 

that: 

^k = ^ y [l(s£V’0 I I Yn)f] {En - Eo) (56) 

n>0 

We shall also use the susceptibility which is given by: 

1 


Ak=2E[l(^£^ol^")l' + l(("k)*V'o|V'n)f 


n>0 


Er,. — Ef) 


(57) 


Clearly Xk < oo due to (541 and the fact that limif„ = d-oo. 
For a positive integer n > 0 let 


dr), .— 


\J\{K'>l^a\Yn)? + l((Sk)*V'0 I V'n)!" 


Y En — Eq 
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bn. ■ — 


l(Sk^O I ^n)|^ + l((Sk)*V'0 I ^n)|^ {En “ Eq) . 


Then = i^anbnf' and (by Schwarz inequality): 


5k = 


(E 


Clribr 


< 




Notice, that since sg is a multiplication by a smooth function and H is a second order differential 
operator, the double commutator appearing in (55) is an operator of multiplication by a smooth 

(58) 


function, so it is bounded. In fact by a direct calculation one gets 




Using this estimate we obtain the inequality: 

5k<Xk-I?k< (59) 

This shows that an upper bound for yt implies the upper bound for g\^. 


3.4 Reflection Positivity arguments 


Let us now proceed along the general line of RP arguments. We perturb the Hamiltonian in the 
analogous manner as it was done in the case of positive temperatures 0, m- To do it, let us first 
modify the original Hamiltonian (43) by a constant (so irrelevant) term introduced by the potential: 


V = ^ ^ [(cos(/3x — cos:/3y)^ + (sintpx — sin(/3y)^] = Vq + J ^ 1 


(60) 


(xy> (xy) 

Now, for a function 6:A9xi—>'6x€C defined on sites, let us consider the perturbed potential V{b): 


(xy) 


and the perturbed hamiltonian 


m = ^ U - cos,,, + (,,t + (sinv, - oinv,)-- 


H{b) :=T + U(6). 


(61) 


(62) 


Let us remark, that if functions b and 5' differ by a constant function then E[{b) = H[b'), in particular 
for b being a constant H{b) = H{0). 

Clearly the perturbed Hamiltonian satisfies all assumptions of Thm|3.ll We will use the Reflection 


(63) 


Positivity and the operator inequality in the form given in (34) to prove the following: 


Theorem 3.5 Let Eg^b) be the ground state energy of the perturbed hamiltonian (62). Then: 

Eo{b) > Eo{0) =: Eq 


A dill): 


The rest of the subsection is devoted to the proof of above theorem. Let us recall the definition of 
A := {y: € : -N + 1 < Xi < N, i = 1,... ,d} (64) 


We divide the system into two identical subsystems Al,Aji so that A = A^ U A^j; Al is an mirror 
image of Aji under reflection in the H (hyper) plane given by the equation xi = |. This way, the subset 
A]^ contains all sites, where the first coordinate is negative or 0, and A^^ ~ sites, where it is positive. 

The Hilbert space TLa of states of the whole system is a tensor product of two spaces: TLa = C^TZ, 
where C (TZ) is a space of states of subsystem defined on Al {Ar). 

For a function b : A —>■ C define: 


( b^ = b on A A, &L is defined on A^ by the mirror symmetry 

\ bn = b on Ar, bR is defined on Ar by the mirror symmetry 


(the standard trick in RP). Our first step to prove inequality (63) is the following : 


(65) 
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Figure ll Division of the system into two identical subsystems Aj,, A/^ by the symmetry plane IT. The illustration concerns open 
boundary conditions. Let us stress that in the paper we have to do with periodic boundary conditions, where analogous division can 
be made. 


Lemma 3.6 Let 6 : A —>■ C &e a function and bL^bn be related to b as in {65). Let hamiltonians 


H{b).,H{bR),H{bL) be defined by [62) and Fo(6),Fo(&i?),F o(6l) denote their ground state energies 
respectively. Then: 

2Eo{b) > Eoibn) + EoibL) (66) 

We move the proof of this lemma to the end of this subsection. 


Now, we will apply RP arguments to show how the inequality ( |66| implies (631. For a given function 
& : A —C let us call a non-zero bond for b a pair of nearest neighbours (x, y) with 6(x) ^ b{y). Clearly, 
if b ^ const then there are I > 0 non-zero bonds for b and we can choose the symmetry plane, which 
crosses at least one of them. Notice that if 1^, Ir are number of non-zero bonds for and b^ 
respectively then, we have: 

II +Ir < 2Z; 

i.e. at least one of the numbers {Il,Ir) is less than 1. Therefore for any non-constant b we can pass 
to a constant one by applying finitely many replacements 6 —> or b ^ bn (for different symmetry 
planes). For a given lattice A the maximal number of steps is bounded independently of 6 - let K 
denotes this bound. 


Our modified hamiltonian (621 is a positive operator, so the set {ifo(&)} is bounded from below and 


let us define E := inf {ifo(6)}. From the inequality ( |M| ) it follows that for e > 0: 

if E < Eolb) < E -f e then E < E^ibif) < E -|- 2e and E < EQ{bn) < E -|- 2e 


(67) 


So let e > 0 be given. There exists b with E < ifo(6) < E -1- 2 ^e; replacing succesively b hy b^ or bn 
and applying (671 we obtain bo = const with 


E < Eo{bo) < E -|- e; 


but since Eo{bo) = Eo{0) — Eq it follows that E = Eg 
prove the lemma 


so to prove the inequality (63) it remains to 


Proof of the lemma 3.6 ' Consider the perturbed hamiltonian El{b) = T -|- V(6), where V{h) is given by 


(62l. The proof consists of two steps: the first is rather involved - by using the lemma (2.81 and prop 
(2.101 we will show that for a smooth function S £ 0 77., there exist ^n, S £ 0 77 such that: 

2 («-1 H{h)-^) > I H{bL)^L) + i^R I Hibn)^R) (68) 


Since, by the variational principle ('I'l | H[bn)'^L) > Eo(^l) and {'i’n \ 77(&fl)'I'i?) > Eo{bn), we get: 

(^-L I H{bL)'fL) > Eo(6l) + Eo{bn). 

For 'i’ = ipo - the ground state of H{b), we obtain the inequality ([M|. 
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To simplify notation let us define a^ib) := costpx ~ b^', we will write just if it is clear what is b. 
Now, V{b) reads: 

= T X! kx - Oyp + (sinv5x - sin(,9y)^ 

<xy> 

Let us define subsets Bl C Al, -B/j C Aj^: 


Bl '■= {x € Al : the first coordinate of x is 0 or — A^ + 1} 
Bji := {x S Afj : the first coordinate of x is 1 or A^} 


(69) 


Notice that Bl (Br) is the subset of those elements in Ar (A^^) which have (some of) their nearest 
neighbours in A^^ (A^). The potential V{b) can be written as: 




I I H" 


(sin^Jx - sin(^y)^ + ^ ^ 


<xy>CL 


2 ^ y I 

<xy>C-R 


+ (sini^x — sin(^y)^+ 


J 


E 


’-^x 


+ (sin(^x - sini^x')^: 


x6Bi 

where x' is the image of x by the reflection across the hyperplane 11. The last term of the sum above 
we write as: 

A E l«x - flx'P + (sin(/Jx - sintpx')^ = ^ E (l«xP + sin^x) + ^ (kx'P + sin^ <Px') + 




xGBl 


:x.eBL 


- ^ (flxflx' + OxOx' + 2 sin (/Jx sin ) 


xeSi 


Let us define functions: 

J 
2 


yL{b)-.= ^ E |ax(&)-ay(&)P + (sin(/9x-sinv3y)2 + ^ ^ |ax(6)P+sin^(/?x 


<xy>CL 


xGBi, 


yR{b) E |ax(&) - ay(6)|^ + (sintpx - sintpy)^ + ^ |ax'(^)l^+sin^ V?x 


<xy>Ci? 


xGBz, 


J 


Vi{b) - ^ E (“xQx' + OxUx' + 2sin(pxSin(px') 


(70) 

(71) 

(72) 


:x.eBL 


Using above notation we can write the potential V{b) as: V{b) = VL(b) + Vh(&) + Vi{b) and the 
corresponding operator V{b) as V{b) = Ul( 6) 0 J + J (g) Vr(6) + Vi{b), where 


Vi{b) := -- ^ (c 


'X 'O *^X 


+ Ox <g) flx' + 2 sin (/9x ® sin yix') • 


(73) 


xGBl 


For the kinetic term T we have T = Tl(^I + I^Tr. 

Now, Let us choose an orthonormal basis T := in C consisting of eigenvectors ofT^: TLi’-y = 
t~tjjj and the corresponding involution Jp. Note that since Tp is self-adjoint it commutes with Jp: 


JpTp Jp = Tp. 


(74) 


Let [/:£—>■ 77. be a unitary operator and (j)^ := Utp^ be the corresponding basis in 77. 

Let us also assume that 

Tii = UTlU* (75) 


That means, in particular, that cj)^ S D{Tji) and Trc/)^ = tycj)^ 


Let 4' =: r(c) G £(g)77 be a smooth function', then it belongs to domains of T*^, (Tp®/)*^, (/(gT^)* 
for k = 1,2,3,... and, using and (|72[), we can write: 


2(4^ I 77(6)4') = 2(4-1 T4') -k 2(4' | U(6))4') = 

= (^ |(Tp®7)^) + (4'|(Tp® J)4')+ (76) 

+ (4'|(/®Tfl)4') + (^|(/®Ti^)4')+ (77) 

-k (4' I (VLib) (g) 7)4') -k (4' I (V^(6) (g) 7)4')-k (78) 

+ (4' I (7 0 t^(6))4') -b (4' I (7 0 t4(&))^)+ (79) 

-b2(4'| V>(6)4') (80) 
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Now, due to the lemma |2.8| and proposition |2.10[ we are going to rewrite various terms appearing in 
this eqnality in different form. 


By formulae (351 and (751 for terms in (76lwe have: 


(4-1 (Tl 0 m) = (f ((7|c|) I (Ti 0 I)tiU\c\) 

(vE-1 (Tz. 0 m = (f (C/|c|) I (/ 0 Tzz)f (C/|c|)) ; 


Using (36l with d = Jtc*Jq, (then |(i| = Jn|c*|Ja and i}{U*\d\) = r(|c*|(7) ) and (75l for terms in 

(0: 


(«-1 (7 0 = (f{\c*\U) I (7 0 Tzz)f(|, 

(«-1 (7 0 Tzz)^') = (t{\c*\U) I (Ti 0 I)f{\c*\U) 

For terms appearing in (781 and (791 by the use of (291, (301, (281 and ( [M| ) we get: 

(vE-1 iVLib) 0 7)vl/) = (f(C7|c|) I (Mb) 0 7)f(7/|c|)) 

I {VUb) 0 7)^) = (f(C7|c|) I (7 0 UMVL{b)rJTU*)f{U\c\)) 
(vP I (7® E^(6))vl/) = (r{\c*\U) I (7® t^(6))f(|c*|C/)) 

(^1(7® Ur( 6))^) = (f (|c*|C/) I {U*Ja{W)rJnU 0 7)f (|c*|C/)) 


Now we are going to use the inequality (|34| to estimate the term in (80). 

Since (5' | (a^ 0 Ox' + Ox <8) aU)^') is real we have: 

('E' I {a^ (g) flx' + 0 x 18 ) oU)'!') < |('I^ I <8) fix' + Ox <8) ^)'E')| < 

< K^-l (a;r(8)ax')^)l + K^-l (ax<8)aU)^')| 


Applying the inequality (34) to each of two terms we get: 

2|(4-|(ax®ax')^)l < (f(C/|c|)|(a,f®C/JraxJrU*)f(U|c|)) 


+ (f(|c*|7/)|(C/Vnax'JnU(8)ax0f(|c*|7/) 


and 


2 I (4-1 (ax (8) oU)^) I < (f(C/|c|)|(ax<8)UJraxJrU*)f(t/|c|)) 

+ (f(|c*|C/)|(C/VnaUJnU®aU)f(|c* 

Similarly, since 2 (4/1 (sintpx <8) sintpx')^) i® 

2 (4^ I (sinvJx <8) sin:/Jx')'I') < 2 |(4' | (sinyix <8 sini^x')'^')l ^ 

< (f(C/|c|)| (sin(^x 8 UJr sin(/9x<7rU*)f ((7|c|)^ + 

+ ^f(|c*|7/) I {U*Jn sm(p:^>JnU 8 sint/Jx^rdc*!?/)^ 

Adding these three inequalities and multiplting by a negative number {—J12) we obtain the estimate 


for (80): 


2(4/1 Vim) > (f (U|c|) I Kl f(C/|c|)) + (f (|c*|U) I Tfflf (|c*|U)) , 
where operators and Kn are defined by: 


J 


Kl ^ [ox 8 C/JrOx JrU* + ttx 8 UJrOx JrU* + 2sin(/3x 8 UJr sin(/3x>/rU*] (81) 


xGBi 


Kii:=-- ^ [U* Jnflx'JnU 8 Ox'+ U* Jnflx'JnU 8 Ox'+ 2U* Jo sinyix'T^nU 8 sin(/5x'] (82) 


xGSi 
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Putting all together we obtain the following inequality : 


2(4-1 > (^L I Hi 4-^) + {'fR I H 2 'fn), where (83) 

^ vl/^;=f(t/|c^, vI/^:=f(|c*|C/) (84) 

Hi :=TL(^I + I(^TR + VL{b)(^I + I®UJr{VLib))*JrU*+K l (85) 

H2 ■.= TL<^I + I<^TR + U*Jn{\^{b))*JnU^I + I^\^{b) + KR ( 86 ) 


We would like to have Hi = H{bR), H 2 = H{bR). Notice that the Hi and H 2 depend only on 
composition UJr (by (27l JqU = UJr)', remember also that we have assumed in (751 that T/j = 
UTlU*. 

Let M : A —)■ A be the reflection across our (hyper)plane 11 and let Um : £ —> 77. be the unitary 
deflned by 


■ ■) := V'(7>xi,7>X2,---) (87) 

y* := M(xj), Xi e Az,, j = 1,..., |Al|, 

and Jo be the complex conjugation (we will use Jq for conjugations on C and 77). It is clear that 

JqUm = UmJq , TlJ o = JqTl , Tfl<7o = JoTzj, 


Define the unitary U := JqUmJt = UmJqJt', then 

UTrU* = c/m Jo JrTi Jr JoI/m = UmTlU*m = and UJr = JqUm = UmJq- (88) 

We will analyze Hi. Clearly we have 1/^(6) = Vz(&l); since the function 1/^(6) is real (Vl(6))* = Vl( 6) 
and JoVl(C))Jo = Pl( 6). Therefore 

U Jr{VL{b))* JrU* = UMVL{b)U*M = VRibr) 


This way we obtain: 


JCi = T + VLibr) ® I + I ® VRlbr) + Kl 


Let us show that Kr given by (811 is equal Viibr) deflned in (|73[): 


C/Jpo-x)^)JpC/ — — C/m(cos:Px bjb)UM — cos:/:^x' ^l(^ ) — 

= ay,i{bL), 

where x' := M(x) (because for x G Az b(x) = &l(x) = 6z(x')); the next term: 


UJrax(b)JrU* = UMaxib)Ulj = ay^i{bL) 


and Anally 


So we get: 


C/Jr sin(/Jx JrC/* = C/m sin (/JxC/m = sin(/?x' 


Kl = -^ ^ |^ax(&L) 0 ax'(6L) + ax(5L) 0 ax'(6L) + 2sin(px C) sin:/Jx 

xSBi 

and this is Vi{br), so really Hi = H{bR); in the same way one gets H 2 = H{bR). 

This way the proof of the inequality (681 and the lemma 3.6 is completed as well as the proof of 
the inequality 


3.5 Estimations giving LRO 


In this subsection we Anally complete the proof of LRO, i.e. we show the inequality (521. The inequality 

d2j;(A6) 


(631 implies that 


dA2 


> 0 


(89) 


A=0 


17 


















for an arbitrary b. It turns out that if we take b being plane wave with the wave vector k then we get 
(52l. Let us present calculations in more details. 


For the moment, let us keep the b function being arbitrary. Write the perturbed Hamiltonian (621 
with the b function rescaled by a factor A G K: 


H{Xb) = H{0) + \H'{b) + 


where 


H'{b) ■= -JRei ^ (cos- cosily)(6x - &y) j , C{b) := ^ ^ l&x - 

\<xy> / <xy> 


(90) 

(91) 


( 2 ) 


Let EoiXb) be the ground state energy of the operator H{Xb), and A Eo{b) - the correction to ground 
state energy in the second order perturbation theory for the Hamiltonian H{0) + XH'{b), i.e. 

p (r\ \{i^n \ H'{b)'tpo)\'^ 

n>0 ^ " 


Therefore: 


Choose now the b function as 


d'^E{Xb) 


dA2 


( 2 ) 

= A Eo{b) + 2C{b) > 0. 


A=0 


bx = 


.,2k-X 


With such a choice (91) reads: 


E = J£(k) 


v1^ 


(92) 

(93) 

(94) 


<xy> 


and, since 


^(cos(^x -cos95y)(e*‘"’' -e*’"'5^) = 4^(2- 2coskj) = 24£(k), 


(xy> 




Therefore by (92): 


H'{b) = -2J£(k)i?e(4), 

-4J^P(k) 5: + 2j£,k) > 0 


n>0 


En — Eq 


By the similar computations, replacing b by ib we obtain: E{'{ib) = 2J£(k)Im(s^) and 

-4j“£yk) y: I^ ^ ^ 


n>0 


En — Eq 


Adding these inequalities we get (remember J£(k) > 0): 

KV'n |i?e(4)V:o)P KV^n I Jto(4)V:o)P ^ 1 

En-Eo En-Eo - J£{k)- 


(95) 


For complex numbers a := {ipn \ SkV'o) and (3 := [ipn \ (Sk)*V'o)) by the paralleogram law : \a + /3p + 
\a - = 2(|ap + |/3p, we get 

KV'n I i?e(s£)V’o)P + KV'n I -fw(4)V’o)P = ^ \ SkV'o)!^ + K^/’n | {st)*'lpo)\‘^) , 

and finally: 

1 y- l(V'n|SkV’o)P , lii’n I (4)*V’o)P < 1 

2^^ En-Eo En-Eo - J£{k)' 
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The LHS of this inequality is just succeptibility Xkj compare (57), so we have Xk < Using (58) 

and (59), we obtain the estimate (52): 


5k < Xk ■ < Xk^ < 


4/ - 4/J£:(k) 


The proof of theorem |3.3| is complete. 


4 Summary 


We have extended the Kennedy-Lieb-Shastry-Schupp matrix inequality to the case where matrices are 
replaced by certain infinite dimensional operators. Similar result has been proven in for another 
matrix inequality - the DLS lemma, which is crucial for the proof of occurrence of LRO in the system 
of interacting rotors in low temperatures in d > 3. 

With the use of this inequality and Reflection Positivity technology, we have formulated sufficient 


condition (44) for ordering in the ground state of the system of interacting rotors. In particular, the 


LRO is present in d > 2 for sufficiently large value of IJ. This way, we have shown the occurrence of the 
LRO in the ground state of interacting rotor systems in a direct manner. In the paper uni, analogous 
result has been proven by Reflection Positivity technique, but without checking some assumptions 
(validity of certain limiting procedure). Our present approach does not suffer from this drawback. 
This result has also been obtained in non-rigorous way by scaling and field-theoretic arguments |12| 

We are convinced that our result can be extended to other rotor systems: other (bipartite) lattices 
and larger space of internal degrees of freedom, for instance, for Oin) systems. 

One can pose the problem concerning the occurrence the ordering in opposite situation, i.e. for 
quantity IJ being small. To our best knowledge, this is an open question. One can suspect that the 
LRO should be absent. Such expectation is motivated by the paper [21], where somewhat similar 
result has been proved: There is no ordering in the anharmonic crystal model provided mass of the 
oscilator is sufficiently small. 

There are numerous interesting rotor-like systems, which do not fulfill conditions allowing an ap¬ 
plication of Reflection Positivity techniques. One of most important of them, is the lattice system of 
interacting bosons (for instance, the Bose Hubbard model). The Hamiltonian of this system, written 
in the language of coherent states, becomes the Hamiltonian of interacting rotors of the form (43), plus 
one term more (see, for instance [55]). This last term spoils the Reflection Positivity, an it seems to be 
not possible to apply these techniques to the analysis of interacting boson systems. (Only exception is 
the paper m, where the Bose-Einstein condensation has been proved for hard-core bosons on optical 
lattice. Here, the term spoiling RP is absent due to the hard-core condition). Here we tackle with the 
long-standing and important problem: How to extend the range of applicability of Reflection Positivity 
technique, which works for certain problems, and does not work for apparently very similar ones. 

Another interesting problem is the occurrence of the Kosterlitz-Thouless transition |53| in the two- 
dimensional rotor system. On physical grounds, one can expect occurrence of this transition, at least 
for large momentum of inertia (the quantum-mechanical rotors should not differ too much from the 
2d XY model, for which such a transition has been rigorously proven |25| 1 However, we are not aware 
on rigorous results for interacting rotor systems. 


Acknowledgments. We thank Przemyslaw Majewski for discussions on early stages of this work. 
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